Scattering of Surface Plasmon Polaritons by one-dimensional inhomogeneities. 
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The scattering of surface plasmons polaritons by a one-dimensional defect of the surface is theoret- 
ically studied, by means of both Rayleigh and modal expansions. The considered defects are either 
relief perturbations or variations in the permittivity of the metal. The dependence of transmission, 
reflection and out-of-plane scattering on parameters defining the defect is presented. We find that 
the radiated energy is forwardly directed (with respect to the surface plasmon propagation) in the 
case of an impedance defect. However, for relief defects, the radiated energy may be directed into 
backward or forward (or both) directions, depending on the defect width. 
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I. INTRODUCTION 

The study of electromagnetic excitations local- 
ized on metal-dielectric interfaces, or surface plas- 
mon polaritons (SPPs), has become of essential im- 
portance due to its potentiality for the implemen- 
tation of sub- wavelength photonic circuits ~£ As a 
result of such study, phenomena like total suppres- 
sion of reflection^ 3 - transformation of polarization^^ 
enhanced light transmission^^ formation of band- 
gap a 9 ' 10 ! 11 ! 12 and other interesting properties of plas- 
monic crystals have been discovered. Lately, a great 
deal of attention has been devoted to the creation of 
optical elements for SPPS) 13 ' 14 ^^^ as well as to the 
efficient coupling of light into and out of SPPs. These 
latter problems require a precise knowledge of the scat- 
tering coefficients of the dispersion centers (i.e. devi- 
ations from a flat metal-dielectric interface) placed on 
the path of SPPs. However, while many theoretical 
works have studied the problem of SPP scattering by 
rough surfaces, scattering from simple geometries is 
not so well known. This was undoubtedly due to the 
previous lack of reproducible experimental data, which 
are now available thanks to recent advances in the con- 
trolled patterning of metal surfaces. 

From the theoretical side, the calculation of elec- 
tromagnetic (EM) fields on a metal surface in the 
optical regime is a well-defined but difficult prob- 
lem. Although the macroscopic Maxwell equations 
govern very accurately the interaction of the EM 
fields with the solid, their solution is difficult due 
to the different ranges of length scales involved 
(system size, wavelength, skin depth etc.). Sev- 
eral techniques have been applied to this problem, 
each of them with their advantages and drawbacks. 
The Greens's Dyadic techniqu o 18 i 19 ' 20 or the Discrete 
Dipolc Approximatio n 21 ' 22 ^ 3 are virtually exact meth- 
ods which suffer from the large (quite often pro- 
hibitive) numerical cost involved with inversion of huge 
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FIG. 1: (Color online) Schematic illustration of the studied 
system: SPP scattering at the inhomogeneity formed by 
either (a) the perturbation of the interface profile h(x) or 
(b) the surface impedance of the metal £(x). 



matrices and the need to calculate cumbersome Som- 
mcrfeld integrals. On the other hand, the mode match- 
ing technique is computationally much simpler, but it 
can only be directly applied to the case of indentations 
in a metal film (and not to protrusions) i 25 ' 26 Addition- 
ally, it requires the use of surface impedance boundary 
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conditions (SIBC), which are only applicable when the 
metallic dielectric constant, e, satisfies |e| ^> 1. An 
alternative approach based on the Rayleigh approxi- 
mation (which is valid for small scatterers) has been 
also extensively applied in the context of EM scatter- 
ing by rough surfaces^ Even within this approxima- 
tion, the calculation of the scattering coefficients re- 
quires solving a difficult integral equation i 28 i 29 i 30 Sev- 
eral works have rendered this integral equation into a 
more manageable form by means of an additional ap- 
proximation, which assigns a geometry-dependent "lo- 
cal impedance" to the surface relief (see, for instance, 
Refs. I3lll32ll33[ ) . Up to our knowledge, all works have 
concentrated on scatterers with translational symme- 
try in one direction, onto which the SPP impinges at 
normal incidence, except for the case of a circularly 
symmetric defect, considered in Ref. [3^ by using the 
reduced Rayleigh equation. 

In this work, we use the Rayleigh approximation 
combined with the SIBC at the metal-vacuum inter- 
face, systematically taking into account the geometry 
of the surface. Our first result is an integral equa- 
tion governing the SPP scattering by arbitrary per- 
turbations, which can be either spatial variations in 
metal permittivity or in surface profile. This integral 
equation is as simple as the one defined by a "local 
impedance" ; 31 i 32 ' 33 but is far more accurate, as shown 
by comparisons with results obtained using the mode- 
matching technique. Due to its simplicity, the method 
developed provides a clear physical description of SPP 
scattering effects. As a first application, here we con- 
centrate on the scattering properties of a single inho- 
mogeneity, as a function of shape and geometrical pa- 
rameters defining the defect. 

II. RAYLEIGH EXPANSION 
APPROXIMATION 

A. Fields representation and boundary 
conditions 

Consider a monochromatic surface plasmon with fre- 
quency cu, impinging along a metal- vacuum interface at 
normal incidence (with in-plane wavevector k p ) onto 
an inhomogencity region, see Fig. [T] In keeping with 
other works on scattering of SPPs, we consider that 
the metal is lossless. From a physicist's point of view, 
we expect this to be valid in the systems we have in 
mind, as far as the dimension of the inhomogeneity in 
the direction of SPP propagation is smaller than the 
absorption length. More mathematically, SPP scatter- 
ing channels arc well defined in a lossless metal, and 
current conservation provides a strong (although not 
definitive!) test for the theory. Nevertheless the pro- 
posed method could, mutatis mutandis, accommodate 
for absorption. In this case, instead of scattering coef- 



ficients, the outcome of the calculation would be spa- 
tial EM field distributions which could be used, for 
instance, to analyze scanning near-field optical micro- 
scope or leakage radiation experiments. 

The considered inhomogencities can be due to vari- 
ations either in the surface impedance or in the sur- 
face profile. Although the integral equations govern- 
ing scattering will have the same functional form for 
both cases, the derivations are slightly different. Let us 
concentrate first in the latter case, which is somewhat 
more involved. Suppose that the metal has a dielectric 
constant e (and correspondingly a surface impedance 
£ = l/-\/e) an d that the surface relief profile has the 
functional form 

z = h(x), h{x) = J dk h(k) exp(ikx). (1) 

Assuming that both the variation of the surface relief 
and its derivative are small (\h\ -C A and \d x h\ <C 
1), we may represent the field over the surface in the 
form of Rayleigh expansion^ The non-zero EM field 
components in the vacuum half-space can be expanded 
in terms of the incident SPP plus scattered field as 

E x (x,z) ( -k zp /g 

E z (x,z) \ = < -k p /g > exp [i(k p x - k zp z)\ (2) 
H y (x,z) J { 1 J 

r f E x (k) ] 
+ dk < E z {k) > exp [i(kx — k z z)] . 
J { H y (k) J 



where g = uijc = 2ir/\ , k z = yj g 2 — k 2 and k zp = 

yjg 2 — k 2 = —£g. The branch of the square root 

should be chosen such that Im(fc 2 ) > 0, in order to 
satisfy the radiation condition. Notice that, within the 
SIBC, the SPP dispersion relation is k p = g(l-£ 2 ) 1/2 , 
which approximates the exact SPP dispersion relation 
k p = g[e/(l + e)} 1 ^ at |e| > I. 

Within the SIBC^ the EM fields should satisfy 

E t (x, z) = £H t (x, z) x n(x), at z = h(x), (3) 

where n = (n x ,0,n z ) is the unitary vector normal to 
the surface (directed into the metal half-space) and 
subscript t corresponds to the tangential components 
of the fields. 

Notice that the SIBC assumes that the radius of 
curvature of the surface is much larger than the skin 
depth. However, as the comparison with results ob- 
tained by the modal expansion method will show, the 
SIBC still represents accurately scattering by defects 
where this condition is not fulfilled at a small num- 
ber of points (as occurs for rectangular and triangular 
indentations). 
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B. The integral equation for the scattered field 
amplitudes 

Expressing the tangential component of the electric 
field as 

E t (x,z) = E(x,z) - n(x)\E(x, z) ■ n(x)], (4) 

the x component of Eq. © gives 

E x (x, z)n z {x)—E z (x, z)n x (x) = £H y (x, z), at z = h(x) 

(5) 

Substituting the fields from Eq. © into Eqs. ([5]), and 
using Maxwell equations E x = —(i/ g){d / dz)H y and 
E z = {t/ g){9 / dx)H y , provides an integral equation for 
the amplitudes of the scattered fields, for instance, for 
the magnetic component, H y (k). A much more man- 
ageable equation can be obtained for smooth surface 
inhomogeneities by expanding the boundary condition 
over the small parametes \d x h\, \h\/\ <C 1. As the sur- 
face normal vector has the form 



n = 



d x h 

yft + (d x hf 



(6) 



the expansion of these vector components up to 
second-order terms gives 



1 



1 - 2^) 2 



o(\d x h\ 4 ), 

d x h + 0(\d x h\ 3 ). 



(7) 



Following an analogous procedure on the exponents, 
exp(— ik zp h), exjp(—ik z h), appearing in the boundary 
conditions (but in the parameter |/i|/A), after some 
straightforward algebra we finally obtain an expanded 
integral equation, in which as many expansion terms 
(in powers of d x h and |/i|/A) as necessary should be 
retained in order to satisfy energy conservation up to 
a required accuracy. However, as will be shown later, 
just considering up to the terms linear in d x h, \h\/\, 
already provides accurate energy conservation. More- 
over, we will show that this even occurs for rectangular 
or triangular defects which are, in principle, unfavor- 
able cases as the shape slopes are not small everywhere. 
From the physical point of view this means that small- 
scale spatial components of the fields do not contribute 
essentially to the scattering. It is convenient to define 
the dimcnsionlcss wave- vector components q — k/g 
(so that q p = k p /g = y/l — £ 2 ) and q z = k z /g, and 
the dimcnsionlcss Fourier amplitude of the relief de- 
fect 77(g) = ig 2 h{k). Additionally, the renormalizcd 
field amplitudes r(q) arc defined by 

r(q)=gH y (k)G(q)-\ G(q) = l/(£ + q z ), (8) 

where G(q) is the Green's function corresponding to 
the unperturbed SPP. 



Then we find that SPP scattering is governed by 
r(q) + f dq' U(q, q')G(q')r(q') = -U(q, q p ). (9) 



In this equation U(q,q') is the scattering potential 
which, in general, can be expressed as a series expan- 
sion in the Fourier image of the defect profile 77(g) as 

U(q,q') = U^q') + U 2 (q,q') + U n (q,q') ~ V n . 

(10) 

Up to first order in \d x h\, \h\/X, the explicit expression 
for the potential is2& 

U(q, q 1 ) = [q'(q - q') - q'Mq'r'Uq - q 1 ). (11) 

We have to keep in mind that, within the "local 
impedance" approximation, the SPP scattering is also 
governed by Eq. © but, in this case, the scattering 
potential is Ui oca i(q,q') = r)(q-q') (e-l)/e ps t](q-q'). 
This different functional form is not irrelevant: as we 
will show later, the fact that U(q p , q p ) = (property 
not shared by Ui oca i) has important consequences in 
the scattering of SPP by surface reliefs. In order to 
ascertain which is the best approximation, subsequent 
sections present comparisons between results obtained 
with the modal expansion and both potentials U and 
Uiocai- Let us anticipate that the comparison favors 
the scattering potential U(q,q') defined by Eq. (fTTj) . 

Another point favoring U(q,q') against Ui oca i(q,q') 
is to consider the perfect conductor limit (e — ► —00), 
where the SIBC must transform to (d/dn)H y (x, z) = 
0, evaluated at z = h(x). In this limit, while Eq. 

transforms correctly, the use of "local impedance" 
leads [by using Eqs. (5) and (25) of Ref. [3l[ to the 
boundary condition (d/dz)H y {x,z) = g 2 h(x)H y (x, z), 
evaluated at z = 0. The root of the problem seems 
to be that the series expansion in h(x) [Eq. (2.23) in 
Ref. [3l| diverges as e — > 00 (note, for instance, that 
the second-order term is inversely proportional to the 
skin-depth, d ~ 1/V^e). 

To complete this section, let us point out that the 
previously outlined formalism can be generalized to 
the case of defects due to impedance inhomogeneities. 
In this case, the metal surface is flat, but the surface 
impedance becomes a function of the x-coordinatc 

£(x) = £ + = £ + J dki(k)exp(ikx). (12) 

After applying SIBC, which in this case reads 
E x (x,z = 0) = £(x)H y (x,z = 0), we find that SPP 
scattering is also controlled by Eq. ([9]). The difference 
is that now the scattering potential is 



U(q,q') = v(q-q'), 



(13) 



where 77(g) = g£(k) is the dimcnsionlcss Fourier ampli- 
tude of the modulation defect. 
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C. The transmission, reflection and out-of-plane 
scattering 

Once the coefficients r(q) are obtained, the integrals 
in Eqs. define the EM field everywhere in vacuum. 
However, in order to find the scattering coefficients, 
only the asymptotic values at large distances from the 
scattering center are needed. It is then convenient 
to prolong the sub-integral functions into the complex 
plane, taking into account the presence of poles. While 
the renormalizcd Fourier image of the field, r(g), is not 
expected to present poles, the Green function G(q) def- 
initely docs. In order to retain its causal character, an 
infimtesimally small damping should be included as 
£ — * ilm(£) + 0. As a result, the magnetic field at the 
surface 2 = may be written in the following form 



H y (x 
H y (x -> 



oo,0) = (1 + r) cxp(ifcpx), 
-oo,0) = exp(ik p x) + pexp(— ik p x), 



where 



) 



(-%»)• 



(14) 



The energy flux scattered out of the metal-vacuum in- 
terface is computed by integrating the Poynting vector 
over a rectangular parallelepiped defined by the plane 
2 = and walls placed in the far-field parallel to the 
planes Ox, Oy, and Oz. Then, taking into account 
that the power per unit length of the incoming SPP 
is g p /(4|£|), energy conservation law has the following- 
form 



1-S-T-R = 0, 



(15) 



where the reflection R, transmission T, and out-of- 
plane S scattering coefficients are 



R=\p\ 2 , T=|l+r| 2 , 



S = 



4^1 



dq-q z \G(q)r(q)\ 2 . 



The total out-of-plane scattered current can then be 
written in the following form 

tt/2 

S = J d6D(6), 

—tt/2 

D {6) = ^ cos 2 9 \r{q(6)]G{q(6)}\ 2 
47r|f|cos 2 6> , , ,...,2 

- abi 1 * 1 - (17) 

where D(6) is the differential reflection coefficient 
(DRC), which provides the angular dependence of the 



radiated energy. In this expression angle 9 is defined 
from the normal, while q{6) = sin 6* and q z (6) = cos 9. 

To summarize this section, the set of equations ©, 
(|16p and (fl~7|) define the scattering coefficients for an 
SPP impinging normally into an arbitrary set of per- 
turbations (with an axis of translational symmetry) in 
a flat metal surface within the Raylcigh + SIBC ap- 
proximations. Such perturbations can be either varia- 
tions in the surface relief (either indentations or pro- 
trusions) or variations in the surface impedance. 



D. Perturbative approach: qualitative 
description of the results 

In some cases, it is useful to estimate the solution 
of Eq.® using a perturbative approach. Taking into 
account the representation of the scattering potential 
in (p~0|) , we seek the solution of Eq. (J9j) in the form of 
a series expansion in 

r(q) = n(q) + r 2 (q) + r n (q) ~ if. (18) 

The first-order Born approximation (FOBA) gives us 

n(q) = -Ui(q,q p ) ~n(q-q p ). (19) 

The second-order term is 

r2(q)=-U 2 (q,qp) + J dq' U 1 (q,q')G(q')U 1 (q' ,q p ). 

(20) 

As usual, FOBA describes the scattering of the inci- 
dent wave with in-plane wavevector, q p , into a plane 
wave with the wavevector q, the momentum difference 
being provided by a single interaction with the defect. 
Therefore the scattering amplitude for this process is 
proportional to rj(q — q p ). This explains the structure 
of the transmission coefficient. As Eq. (fT6]) shows, T 
has contributions from both the incident SPP (with 
unit amplitude) and from the scattering of the inci- 
dent plasmon "into itself" [the term proportional to 
77(0)] . Analogously, in FOBA the amplitude of the re- 
flection is proportional to rj(—2q p ) and represents the 
scattering of the incoming SPP into the backwardly- 
propagating one (having the wavevector — q p ). 



III. SPP SCATTERING BY A SINGLE 
DEFECT. 

In this section we consider the dependence of scatter- 
ing properties by single ID defects of different shapes. 
The defect width is a and the defect depth is w (see 
Fig. 2). Notice that, within the chosen coordinate sys- 
tem, indentations are characterized by w > 0, while 
protrusions have w < 0. The calculation is performed 
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by solving numerically the integral equation af- 
ter applying an appropriate discretization in q-space. 
As the results are in good agreement with those ob- 
tained within first-order Born approximation, analytic 
expressions for the scattering coefficients are provided 
in several cases. 



A. Scattering coefficients 

Figure [3 renders the results of numerical calcula- 
tions for R, T, S in different instances, fn (a) the 
comparison of these quantities for single indentations 
of different shapes is presented. In the case of a rect- 
angular shape, R behaves periodically with respect to 
a I A, while in the case of a gaussian shape R possesses 
only one maximum. Analogously, the transmittance 
presents oscillations in a/ A, in contrast with the single 
maximum that appears for the case of gaussian shape. 
In the same figure we present also the calculations for 
the triangular shape. These results show that even 
shallow defects present very rich shape and spectral 
dependences. 

Figure [U(b) shows the scattering coefficients for 
both indentations and protrusions of gaussian shape. 
As can be seen from this figure, the scattering prop- 
erties of indentations and protrusions are very similar 
for shallow defects. Nevertheless, for the considered 
parameters, protrusions present a slightly larger cross 
section, resulting in larger values for both R and S and 
smaller ones for T. We stress that, although the inte- 
gral equations where derived assuming that the surface 
was smooth, energy conservation was fulfilled with an 
accuracy better than 1% of the minimum value of R, 
T and S, even for sharp defects. Additionally, we have 
checked that retaining second-order terms in 77(g) in 
the scattering potential leaves the results virtually un- 
altered. Moreover, solving the integral equation with 
the "exact" (not expanded) right-hand side does not 
produce any significant variation in the calculated scat- 
tering coefficients. 

In order to further validate the above mentioned ap- 
proximations, leading to the results presented in Fig. [2J 
additional calculations were carried out with the modal 
expansion technique^. While also using SIBC, this 
technique is applicable for indentations of any depth, 
going beyond the Rayleigh expansion. The compari- 
son is presented in Fig. [5] (a), for the case of a rect- 
angular indentation. As can be seen, the agreement 
is very good, the difference being attributed mostly 
to the fact that, in the modal expansion, ideal metal 
boundary conditions were used for the vertical "walls" 
of the rectangular indentation. We note in passing that 
using the "local impedance" scattering potential gives 
very different results: the values of T and S differ by 
more than an order of magnitude from the ones pre- 



sented in Figs. [2] and [3] (we performed the calculations 
with the potential presented in Ref. [H for the same 
set of parameters). 

We find that FOBA provides an accurate descrip- 
tion for the behavior of the calculated scattering co- 
efficients. This is illustrated in Fig. [5] (b), where the 
FOBA results for the gaussian defect are compared 
with the full solution of Eq. @. Notice that FOBA 
predicts the same reflectance for both protrusions and 
indentations, since R ~ w 2 . Further analysis shows 
that the second order approximation already accounts 
for the small differences between the scattering coeffi- 
cients of indentations and protrusions found in Fig. [2] 
(b). This occurs because, in second-order approxima- 
tion R ~ |w + w 2 • ip(a)\ 2 , where the complex function 
ip{a) depends upon both shape and longitudinal size 
of the defect. 

Besides, FOBA enables us to find analytic expres- 
sions for the reflectance of SPPs by a single defect quite 
easily, in terms of |r/(— 2q p )\ 2 . For instance, for a rect- 
angular or a gaussian defect of width a, we obtain 



■nil) 



Rect 



. A sin(qfi) 
lA— ^-A 77(a) 



Gauss 



iA 



nq lypH 

(21) 

where a = air/X and A = gw for relief defects (for 
impedance defects A corresponds to the maximum 
value of |£|). Therefore, the reflectance of a single relief 
defect can be expressed as 



R 



Gauss 



47r|£| 2 g 2 A 2 & 2 e- 2 ^ a2 



(22) 



R Rect = i 6 |£| 2 2 A 2 a 2 sinc 2 (2^6) . 



Thus, for rectangular defects R behaves period- 
ically as a function of a/ A, possessing minima at 
a/A = n/2q p , n = 1,2,...; while the reflectance for 
defects of gaussian shape presents only one maximum 
at a/A = \/\/2nq v . This is in an excellent accordance 
with the strict numerical solution of Eq. ©, see Fig. [3] 
(a), (b). However, as the transmittance in FOBA does 
not depend upon the width of the defect, higher or- 
der terms in the Born series are required in order to 
reproduce this dependence appropriately. 



B. Out-of-plane radiation due to a single defect. 

In this section we analyze the angular distribution 
of the the energy radiated out of the plane after scat- 
tering. This behavior is represented in Fig. [3J which 
renders the radiation diagrams for SPP scattering by 
gaussian indentations [panel (a)] and rectangular in- 
dentations [panel (b)] , for different defect widths. The 
surface impedance is that of silver at A = 600nm 
(£ = — 0.277i). The case of an impedance defect 
(with zero impedance, as for perfect conductors) is also 
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FIG. 2: (Color online) The dependency of the transmittance, T, reflectance, R, and emittance, S, upon the dimensional 
width of the defect of different shapes in the silver surface (at 600nm). The amplitude is w/A = 0.02. In (a) solid, dashed 
and dash-dotted curves correspond to the rectangular, gaussian and triangular shapes of indentation. Round markers 
correspond to the calculations using mode expansion. In (b) dashed (solid) curve corresponds to a gaussian indentation 
(protrusion); by the rectangular markers the Born approximation is presented. 



shown in Fig. [3] (c). For comparison, the calculations 
for rectangular indentations were also performed with 
the modal expansion method [circles FigOJb)]. Again, 
the agreement between the two methods is quite re- 
markable. Fig. [3] shows that the radiation diagrams 
present a non-trivial dependence on defect shape. For 
impedance defects, the angle at which maximum out- 
of-planc radiation occurs always points along the di- 
rection of propagation of the incident SPP. However, 
for narrow gaussian relief defects (a < A/2), the max- 
imum in DRC occurs at negative angles. In this case, 
as a increases, the angle of maximal radiation shifts 
from negative angles to positive ones. This behavior 
of gaussian defects has also been reported for surface 
plasmons in a thin metal film, excited in a Kretschm- 
man configuration^!. For rectangular defects, the DRC 
behavior is more complex. At small a/A the DRC 



presents one emission lobe at a negative angle; as a 
increases, the emission lobe moves to the normal di- 
rection and after a transition point (at approximately 
a/A = 1/2) the main lobe moves to positive angles, 
while a second emission lobe appears at a negative an- 
gle. Finally, for a/Awl the amplitudes of the two 
lobes are comparable, see Fig. [31 

Within FOBA, the DRC can be analytically com- 
puted. From Eq. (fl~9)) and Eq. (fl~7|) we obtain 

n iq\ 47r %>l£l( shl 9 - 1p) 2 cos2 d \ i - a \|2 
SiM = | e | 2 + C0S 2 g M^9-a p )\ . 

(23) 

As illustrated in Fig. [3] (a) for the case of gaussian 
indentations, FOBA results (square symbols) provide 
and excellent approximation to the full solution. This 
allows us to connect the DRC with the potentials cor- 
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FIG. 3: (Color online) DRC, D(0), for SPP scattering along 
a silver surface (A = 600nm), for different defect types and 
widths, (a) gaussian indentation, (b) rectangular indenta- 
tion and (c) impedance step (the metal is considered ideal 
inside the defect). Dash-double-dotted (black), dashed 
(red), dotted (green), dash-dotted (blue), solid (magenta) 
correspond to defect widths a/A = 0.1, 0.25, 0.5, 0.6, 0.8, 
respectively. The depth of the indentation is w/A = 0.02. 
Squares in (a) correspond to calculations within the first 
order Born approximation. Round markers in (b) corre- 
spond to the calculations using modal expansion 



responding to the SPP scattering within such approx- 
imation. 

For the case of impedance defects, Eq. (jX3J) and 
Eq. @ imply that the renormalized spectra of the field 
r(g) is proportional to rj(q — q p ), which is the Fourier 



image of the defect shifted by the SPP wave-vector—. 
Conversely, for relief defects, r(q) is proportional not 
only to r](q—q p ) but also to q — q p . This is illustrated in 
Fig-Hl for gaussian (left panels) and rectangular (right 
panels) at a/A = 0.8 (top panels) and a/A = 0.1 (bot- 
tom panels). In all panels the modulus of the shifted 
Fourier image of the defect \rj(q — q p )\ is represented 
by the green dashed lines. Fig. Q] also shows the cor- 
responding |r(g)| for relief defects, obtained by solv- 
ing the integral equation @ either strictly or within 
FOBA approximation. 

As expected, the narrower the defect, the more ex- 
tended |7y(g — q p )\ (and therefore |r(g)|) in g-space. 
It is also clear that \r)(q — q p )\ is a smooth function 
for gaussian defects, while it presents oscillations for 
rectangular ones. These two properties are extremely 
useful in order to understand the behavior of DRC. 
Notice that D{9) essentially follows the dependence 
of |r(g)| 2 (see Eq. [T7] within the range — 1 < q < 1 
i.e., the region corresponding to radiative modes, rep- 
resented by the blue shaded areas in Fig. [4J. For 
impedance defects |?/(g — q p ) \ presents its maximum at 
g > 0, so the maximum out of plane emission always 
occurs in the forward direction, exactly in the same 
way as in the model of the local surface impedance 
approximatio n 31 ' 32 ! 33 . However, for relief defects, for- 
ward direction emission is further inhibited by the 
presence of the "cutting function" q — q p in the ex- 
pression for r(g). For narrow defects this results in 
|r(g)|'s having more weight on the g < region, so the 
emission is directed backwards. For wide gaussian de- 
fects (a > 0.5A) the width in q-space of |??(g)| is smaller 
than q p (recall that q p « 1) and the emission occurs in 
the forward direction. For rectangular defects, as the 
defect widens, more and more oscillations in \r)(q — q p )\ 
enters into the — 1 < q < 1 region, leading to different 
emission lobes, as shown in Fig. [3Jb). 

The angles of maximum out-of-plane emission can 
be obtained in FOBA by solving d e Di{6) = 0. While 
in the case of a rectangular shape this equation is tran- 
scendent with respect to variables 9 and a/ A, it may 
be solved explicitly for gaussian shapes. Thus, we find 
that the maxima condition for gaussian reliefs is 



7r 2 (sin( 



where 



<p(9) = 



2|£| 2 tan6» 



(|e| 2 + cos 2 6») (g p -sin6») cos6>' 



(24) 



(25) 



while in the case of an impedance defect this condi- 
tion is much simpler, a/X{9) = ^Jtp(9). These depen- 
dencies are presented in Fig. [5l reflecting that while 
impedance defects emit preferentially in the forward 
direction, emission maxima occurs in the backward di- 
rection up to a I A « 0.42 for the relief case . 
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Dimensionless wavenumber, q 



FIG. 4: (Color online) Dependence of the modulus of the renormalized Fourier image of the scattered field, |r(q)|, upon 
dimensionless wave-vector, q. Left panels show the calculation for gaussian indentations [a/A = 0.8 in panel (a), a/A = 0.1 
in panel (b)], while right panels correspond to rectangular indentations [a/A = 0.8 in panel (c), a/A = 0.1 in panel (d)]. 
The depth of the indentation is w/A = 0.02. In each panel the squares represent the \r(q)\ calculated within first order 
Born approximation, whereas dash-double-dotted (black) and solid (magenta) lines stands for the strict numerical solution; 
the green dashed curve renders the modulus of the Op-shifted Fourier transformation of the defect, \jj(q — q P )\. Notice the 
dips of all curves at g = q p due to the presence in the potential of the "cutting" function \q — q p \. The blue shaded areas 
denote the regions in q corresponding to out-of-plane radiative modes. 



IV. CONCLUSIONS 

We have developed an approximate method (based 
on Raylcigh expansion and surface impedance bound- 
ary conditions) for the study of SPP scattering at ar- 
bitrary ID inhomogeneities, and applied it to the case 
of single defects. 

We have compared the numerical solution of the in- 
tegral equation in fc-space with the calculations made 
by using mode expansion technique. The excellent 
agreement between them (together with the fulfillment 



of the energy conservation law) indicates that the theo- 
retical formulation is correct and much more accurate 
than previous approximations based on local surface 
impedances. 

The case of a single scattering center is analyzed in 
detail. We have compared the scattering by impedance 
and relief inhomogeneities of different shapes, consider- 
ing both protrusions and indentations. We have shown 
that the transmission, reflection and out-of-plane scat- 
tering, are defined essentially by the spectral proper- 
ties of the inhomogeneity. We have shown that, out- 
of-plane radiation after scattering by impedance inho- 
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width /wavelength, afk 



FIG. 5: (Color online) Dependence of angle of maximum 
out of plane radiation as a function of a/X for gaussian 
defects. The case of relief defects is represented in the 
red continuous curve, the blue broken curve corresponds to 
impedance defects. 



mogcncitics is always directed in the forward direction 
(with respect to the SPP propagation). On the con- 
trary, in the case of relief defects, the radiated energy 
may be directed into backward or forward directions 
(or both, in the case of rectangular defects), depending 
on the defect width and shape. 

Further theoretical work will be aimed at the scatter- 
ing of SPP by the multiple scatterers, and, also, at the 
generalization of this approach to 2D inhomogeneities. 
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